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Oh' The stationary one-dimensional tight-binding Schrodinger equation with a weak diagonal long- 

range correlated disorder in the potential is studied. An algorithm for constructing the discrete 
binary on-site potential exhibiting a hybrid spectrum with three different spectral components (ab- 
solutely continuous, singular continuous and point) ordered in any predefined manner in the region 
of energy and/or wave number is presented. A new approach to generating a binary sequence with 
the long-range memory based on a concept of additive Markov chains (Phys. Rev. E 68, 061107 
(2003)) is used. 
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PACS numbers: 72.15.Rn, 03.65.Bz, 72.10.Bg 



In recent years, the problem of transport of electro-magnetic waves (and various excitations in solids) in one- 
dimensional (ID) systems with the long-range correlated disorder has attracted much attention dHQ. The important 
significance of this problem is due to exciting results that revise a commonly accepted belief that any randomness (no 
matter how weak the randomness is) in ID structures results in the Anderson localization £|. R was also believed 
that ID systems could not display a complex dynamic feature such as the metal-insulator transition, which gives rise 
to an appearance of the mobility edges. However, in Refs. [HQ, a highly nontrivial role of correlations was shown, the 
divergency of the localization length for some specific values of energy was observed. Using a perturbative approach, 
it was also shown in Ref. ,1] that the position of the mobility edges and the windows of transparency can be controlled 
by the form of the binary correlator of a scattering potential, which was supposed to be continuous, of a long-range 
and Gaussian. If the correlations are of a long-range, i.e. they decay by the power law, a continuum of extended states 
£^ , may appear in the energy spectrum. In other words, the long-range correlations are necessary for the observation of 
• the metal-insulator transition and co-existence of different type of spectra. A quite simple method was used therein 
£f) \ to construct a random potential with a long-range correlation. From the experimental point of view, the importance 
t — ■ of recently obtained results may be explained by a strong impact upon the creation of a new class of electron devices 
and electromagnetic waveguides, which can be used as window filters having new transmission properties. 

In this paper we point out a new method for constructing the long range correlated sequences of a two- valued site 
potential e(n) with a given correlator and prescribed probability distribution function (PDF), not only Gaussian. For 
this potential the method given in Q does not work. An attempt to construct a correlated dichotomous sequence 
where the metal-insulator transition can be observed was made in Ref. Q. Later, Ref. [U, this result was retracted. 
We study a stationary one-dimensional tight-binding Schrodinger equation (the Anderson model 3), 
■ 

t/'n+i + tpn-i + (E - e(n))ip n = 0, (1) 

o : 

O ■ with a site potential e(n) taking on two different values £o an d £j.. 

Equation is a prototype model describing a propagation of excitations (electromagnetic waves, electrons, photons, 
phonons) in deterministic ordered or random disordered systems (solids or layered super-lattices). The type of 
ordering, i.e., the correlations in the site potential e(n) determine the spectrum of excitation. The wave functions 
of excitations in such systems are usually characterized by the Lyapunov exponent A. Some typical examples of 
sequences having different types of spectra are (spectrum is a set of allowed energies of excitations in the system) : 

1. A system with a regular periodic variation of the site potential e(n) possess an absolutely continuous (AC) 
spectrum and describes extended, delocalized states having bounded wave functions with the Lyapunov exponent 
A = 0. 

2. A system with a random non-correlated potential e(n) has the allowed states with the exponentially localized wave 
functions and displays a point (discrete) spectrum with the positive Lyapunov exponent (the Anderson localization) . 
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3. A potential constructed using the deterministic quasi-periodic Fibonacci chain (see the definition given below, 
Eq. (TUj) exhibits a singular continuous (SC) Cantor-like spectrum characterized by the power-law localized wave 
functions with A = 0. The singular continuous spectrum differs essentially from the absolutely continuous one. This 
spectrum corresponds to a singular continuous integrated density of states (the number of states having the energies 
smaller then a given one). The integrated density of states is the continuous function of energy, but its derivation 
equals to zero almost everywhere. In terms of the forbidden and allowed energy zones, this means that the spectrum 
consists almost completely of forbidden gaps with the infinite number of allowed energies between them. In other 
words, the spectrum of quasi-periodic potential sequence possesses a fractal structure. 

In each of these structures only one pure type of spectrum is presented. A nontrivial result of co-existence of 
continuous and point spectra in one dimensional chain with correlated disorder was presented in Ref. [lj. 

Correlated properties in the site potential e(n) are determined by the pair correlation function, 



1 N 

C{r) =< e{n)e(n + r) > s - < e(n) < /(e(n)) > s = Jim — V /(e(n)), (2) 

N^oo iv — * 



n=l 



where the Cezaro average < • > s = limjv— >oo j? ' can ^> e treated as spatial one. In the Born approximation 

the Lyapunov exponent is expressed, see 1 1 1 1] - in terms of the Fourier transform K of this two-point 
correlation function, 

(e 1 -e ) 2 K(2k) C(r) 

A ( E > = on ■ 2, > k ^ -T.T . K ( r > = TvnV ( 3 ) 

32 sin k C(UJ 
The correlation function K(r) and its Fourier transform are connected by the following relations, 

K(k) = 1 + 2 V KM cos(fcr), KM = - K(k) cos(kr)dk (4) 

Here the evenness of the functions K(r) and K(k) is taking into account. The energy E of the eigenstate and the 
wave number k in Eq. in zero approximation on the strength of disorder arc related by the simple formula, 

£ = -2cosfc. (5) 

Thus, there exists the relation between the correlations in the on-site potential K(r) and localization properties of 
eigenstates expressed in terms of the Lyapunov exponent A(E) (or by means of integrated density of states). This 
observation enables one to construct random correlated sequences with prescribed spectral properties and provides a 
recipe for designing filters of arbitrary complexity. Starting from a desirable spectral dependence A(_E) (or integrated 
density of states) we have to solve the inverse problem of construction a sequence of "symbols" e(n). This program 
was partially implemented in Ref. , where the one-dimensional tight-binding Schrodinger equation with two kind of 
spectra (absolutely continuous and point) and Gauss distribution of the site potential was examined. However, it is 
known [12| , that in the general case the space of states can be decomposed into the direct sum of subspaces with the 
wave functions belonging to the point, singular continuous and absolutely continuous parts of spectra. This statement 
is closely connected with a possibility for presenting an integrated density of states, as for any monotone function, 
in the sum of three functions: the stepwise, singular and absolutely continuous ones. Here we point to a general 
method, which is the extension of that used in Ref. []], and construct a particular example of ID system having all 
three above-mentioned types of spectra [lflj . 

To this end we use a new instrument of constructing the correlated sequence of s(n) with a given correlation function 
K{r) using an additive Markov chain |!4lll5|. Here we give a short description of this method. 

Consider a homogeneous binary sequence of symbols, e(i) = {eq, ei}, ieZ = —2, —1, 0, 1, 2, .... To determine the 
N-step Markov chain we have to introduce the conditional probability function P(e(i) \ s(i — A), s(i — A + 1), . . . , s(i — 
1)). It is a probability of occurring the definite symbol (for example, e(i) — Ei) after A-word Tj^i, where Tjy,, 
stands for the sequence of symbols e(i — N), e(i — N + 10), . . . , e(i — 1). The additive Markov chain is characterized 
by the conditional probability function of the form 

N 

P{e(i) = e 1 | T N<i ) =Pi+J2 F ( r Mi -r)-<e >). (6) 

r=l 

We refer to the amount F(r) (introduced first in Ref. 0|) as the memory function. It describes the strength of 
influence of previous symbol e(i — r) (r = 1,...,N) upon a generated one, s(i). Here p\ is the relative part of 
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symbols £1 among the total number of symbols in the whole sequence (pi is not necessarily equal to 1/2). There is a 
single- valued relation between the memory function F(r) and the correlation function K{r) of the Markov chain, see 
Ref. dHH, 

N 

K{r) = {£ 1 -e )Y d F (r')K{r-r% r > 1. (7) 

r' = l 

This equation allows finding the memory function F(r) and effective constructing the Markov chain with the obtained 
conditional probability P(. | .) introduced in Eq. JBJ. 

Let us present a " Recipe" and enumerate the consecutive steps of the general method for constructing the correlated 
sequence of e n with desired spectra: 

1. Dividing the interval of energy E (or, that is the same, of the wave number k) to the regions where we want to 
have different types of spectra. 

2. Prescribing A or/and K(2k) in these intervals, we construct "by hand" the total Fourier transform of correlation 
function K(2k) for all values of k 6 [— it, ir]. 

3. Calculating K(r) as inverse transform of K(2k). 

4. Solving Eq. 10, we determine the memory function F(r) and the conditional probability function P(. | .). 

5. Constructing sequence of Eq and £\ with obtained function P(. | .). 
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FIG. 1: The line y = ipx + £ "generates" the Fibonacci sequence (0. The strip bounded by two lines y = ipx and y = ipx + 1 
corresponds to an ensemble of all possible Fibonacci sequences with < £ < 1. The line y = ipx + £ intersects a horizontal line 
y = 1 between and 1 (and so does another line y between r and r + 1) at point a (and b) and gives (0) = 1 (and $^(r) = 1). 
Such values of £ give a nonzero contribution to the integral 1121 . 

To construct the total correlation function K{r) we have to calculate the "patrial" ones. Let us present here some 
simple examples of the correlation functions corresponding to the above-mentioned sequences exhibiting the pure 
spectra of periodic, random and quasi-periodic chains. 

It follows from definition @ that the correlation function of a chain of alternating potentials e(n) = ...eqEiSoEi... 
has the form: 

CaAr)= { -^(ei-So) 2 . (8) 

It follows from the same definition (J2J that the correlation function of the sequence with the same potentials 
e(n) — ...EqEoEo... is equal to zero for all distances r. Nevertheless, the ratio C(r)/C(0) — K(r) should be defined as 
unity. Only in this case Eq. gives the correct value of the Lyapunov exponent. 

The correlation function of the non-biased non-correlated random sequence of s n is 

C p (r) = ( £l - £o ) 2 | J/ 4 ;^ = ° (9) 
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Let us derive the correlation function of the quasi-periodic chain of potentials e„ constructed using the Fibonacci 
chain, 

<^(r) = [(r + lV + £]-N + £]. (10) 

as. 

£r(£) = eo + (£i-£o)$e(r), £e[0,l), (11) 

Here ip is the so-called golden number, ip = (V5— l)/2, and £ selects some concrete sequence from among all possible 




r 

FIG. 2: The function K s (r) — C s (r)/C s (0) determined by Eq. 113H vs. the continuous variable r for eo = 0, £1 = 1. 

Fibonacci sequences. The geometrical meaning of this construction is explained in Fig.^ To determine the correlation 
function of this sequence it is convenient to use an ergodic theorem |18| about the homogeneity of distribution of 
y n = n<p — [nip] on the interval < y < 1 for any irrational number ip and calculate correlation function C (r) using 
ensemble average < • > a instead of space average 

C(r) =< e n e n+r > a - < e„ >l < /(£„) > a = / f(e n (0)d£- (12) 

Jo 

From a simple geometric consideration, see explanation in Fig.^ we obtain the following result, 

!2ip-l - {tpr}, {<pr} <l-(p, 
3p-2, l-p<{pr}<p, (13) 

2ip-2 + {ipr}, p< {tpr}. 

Here {x} is the fractional part of a;. The function K s {r) = C s (r)/C s (0) of continuous argument r is shown in Fig. El 
Actually, the correlation function is defined for the integer arguments, which values are incommensurable with the 
period 1/ip of function C s (r) . The function C s (r) of the discrete argument is shown in Fig. [21 Numerical calculations 
of C s {r) carried out using spatial averaging © exhibit a good agreement with the analytical one Below we 

construct an example of the sequence with a prescribed spectral property. Let us suppose that the Schrodinger Eq. 
has all three kind of spectra. Let kn and ka be the wave vectors corresponding to the top and bottom edges of the 
bands where the singularly continuous, i = s, absolutely continuous, i = a, and point (discrete), i = p, spectra take 
place. In fact, we do not need to calculate Fourier transforms K(k) of the correlation function as it is indicated above 
in the paragraph 3 of the " Recipe" . We can calculate the Fourier transform of the composed correlation function 
K(r) presented in fc-space via characteristic functions X[fcii,fc i2 ] V s ) °f * ne intervals [ku, fc^], 



K(k)= E Xlk^m^ik), 

i—{a,s,p} 



(14) 
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FIG. 3: The correlation function K a {r) = C s {r)/C a {Q) of the Fibonacci chain lllfl for eo = 0, ei = 1. 



where characteristic function X[ki,k 2 ] CO is equal to 1 if k G [hi, and equals otherwise. Straightforward calculations 
yield: 



K{r) = 



k t 2 - hi 



Kj(r) + - > iij(r + r ) + iQ(r - r ) ; 



(15) 



Here the functions (sinr'fci2 — shxr'kn)/r' are the Fourier transforms of the characteristic functions. In the case i = p, 
the contribution of the point spectrum to Eq. (|15fl is reduced to one term only, 



sin rk P 2 — sin rk v \ 
k k nr 



r = 0. 



(16) 



Note that in deriving Eq. I|15|) we did not use the concrete form of correlation functions Eqs. JSJ, © and (|13fl . In 
Fig-EJthe result of numerical simulation of the normalized Lyapunov exponent A'(i?) = 32 sin 2 (fc) K{E)/{e\ — £o) 2 is 
presented for the system exhibiting two types of spectrum, namely, absolutely continuous with A' = 0at0<£'<l 
and point spectrum A' = 1.5 at 1 < E < 2 [l9| . 
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FIG. 4: The "normalized" Lyapunov exponent A' of the binary correlated system vs. the energy E. The solid line is the 
prescribed characteristic of system with the stepwise Lyapunov exponent A'. The dots are the result from the calculation of the 
Lyapunov exponent of the sequence e(n) constructed by means of the memory function F(r) obtained by numerical solution 
of Eq. Q for the correlation function K(r) = sin(7rr/3)/(7rr/3), corresponding to the stepwise Lyapunov exponent. 



Thus, we have proposed a method for constructing ID sequences of sites e(n) with a correlated disorder exhibiting 
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a hybrid spectrum with three different spectral components ordered in any predefined manner in the region of energy 
E and/ or wave number k. 

Authors are grateful to L. A. Pastur for the discussion which was helpful in the formulation of the problem under 
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